function dudv = Regulation_Function( t, uv)
dudv=zeros(4,1);
dudv(1) = -uv(1)+uv(3)/(1+uv(2)^uv(4));
dudv(2) = -uv(2)+uv(3)/(1+uv(1)^uv(4));
end

